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Generalizations of Onsager reciprocity relations are established for the nonlinear response coeffi- 
cients of ballistic transport in the effusion of gaseous mixtures. These generalizations, which have 
been established on the basis of the fluctuation theorem for the currents, are here considered for 
mass separation by effusion. In this kinetic process, the mean values of the currents depend non- 
linearly on the affinities or thermodynamic forces controlling the nonequilibrium constraints. These 
nonlinear transport effects are shown to play an important role in the process of mass separation. 
In particular, the entropy efflciency turns out to be significantly larger than it would be the case if 
the currents were supposed to depend linearly on the affinities. 

I. INTRODUCTION 

For long, the understanding of nonequilibrium processes has been limited to the linear regime close to the ther- 
modynamic equilibrium with results such as the Onsager reciprocity relations [1]. Nowadays, great advances have 
been carried out since the discovery of time-reversal symmetry relations for large-fluctuation properties on the basis 
of dynamical systems theory [2H1]- Recently, these relations have been extended to stochastic processes, as well as to 
open quantum systems [5HS] . One of the most detailed versions of such relations is the fluctuation theorem for all the 
currents flowing through a nonequilibrium system [9l 110) . Thanks to this fluctuation theorem, relationships have been 
obtained that generalize the Onsager reciprocity relations from linear to nonlinear response coefficients [llH13j . These 
relations have already been applied to linear chemical reactions , to Brownian sieves and molecular machines |15j , 
and to electronic quantum transport [TBI HZI • 

The fact is that many nanosystcms typically develop highly nonlinear responses to multiple nonequilibrium driving 
forces, as it is the case in electronic transistors for instance |12j . In such devices, several currents are coupled together 
so that different driving forces may induce the response of one particular current. Since they only hold in linear 
regimes, the Onsager reciprocity relations are of limited use to study such nonlinear devices and we may expect that 
their recently discovered generalizations to nonlinear regimes would be valuable to understand the coupling between 
transport properties in far-from-equilibrium systems |11H13) . 

The purpose of the present paper is to address this issue for the process of mass separation by effusion. In this 
process, a mixture of two or more gases is forced to flow from one reservoir to another through a small orifice with a 
size smaller than the mean free path ^18, . Effusion couples the heat and particle flows and kinetic theory has shown 
how these flows depend on the mass of the particles [T5H^ . For a binary mixture, three currents are thus coupled 
together: the heat flow, the total mass flow, and the differential flow of one species with respect to the other. Several 
coupling effects are possible between these three flows and, in particular, the separation between particles of different 
masses induced by either the temperature or the pressure difference between both reservoirs. 

Although a fluctuation theorem has recently been obtained for effusion p2^24j . neither its extension to mass 
separation nor its consequences on the nonlinear response coefficients have been reported on. These coefficients are 
worth to be studied because they obey remarkable relationships that are the consequence of microreversibility [11|- 
113] ■ as the Onsager reciprocity relations are in linear regimes [Ij. Furthermore, these coefficients control the coupling 
between the different currents far from equilibrium and they influence the efficiency of mass separation. In this 
context, a concept of entropy efficiency was introduced as the ratio between the contribution to entropy production 
caused by mixing of the two constituents over the one caused by heat flow |25[ I26j . It turns out that the entropy 
efficiency may be significantly higher in nonlinear regimes than in linear ones, as we show here below. 

The paper is organized as follows. In Section|llj we present the fluctuation theorem for all the currents in the effusion 
of gas mixtures. The generating function of the counting statistics of the different flows is obtained from kinetic theory. 
In Section [lll[ the mean currents and their statistical cumulants are calculated from the generating function and we 
verify that the nonlinear response coefficients indeed obey relationships previously reported in Refs. |llH13j and which 
find their origin in microreversibility. In Section |IV[ the results are applied to the mass separation of two species. 
The entropy efficiency is shown to be significantly higher than possible if linearity was assumed. The conclusions are 
drawn in Section fVl 
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II. FLUCTUATION THEOREM FOR CURRENTS IN THE EFFUSION OF GAS MIXTURES 



A. Thermodynamic formulation 

We consider two gaseous mixtures at different temperatures, pressures and concentrations in two reservoirs TZ and 
TZ' separated by a wail, through which a small pore allows the flow of particles. The reservoirs are assumed to be so 
large that the system is maintained out of equilibrium in a steady state. Alternatively, the far ends of both reservoirs 
can be closed by pistons moving in such a way that the thermodynamic conditions of the gaseous mixtures are kept 
stationary (see Fig. [T]). 



T, P, {«,} 
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FIG. 1: Schematic representation of the stationary effusion of gaseous mixtures through a small pore between two reservoirs 
TZ and TZ' where the thermodynamic conditions are fcept fixed by moving the pistons. These conditions are defined by the 
temperature T, the pressure P, and the particle densities {nij^^i. V, E, and 5* denote the volume, the energy, and the entropy. 
The prime refers to the same quantities in the right reservoir. 

The gases are supposed to be ideal mixtures of c monoatomic species so that the particle densities {n^}^^]^ and 
{n-i}i=i are related to the pressures and temperatures by the perfect-gas equation of state: 



P ^^n,kT 

1=1 



and 



where k is Boltzmann's constant. 

The currents of energy and particles from the left reservoir to the right one are defined by 



Je = -E = 
J,. = -N,. = 



-N' 



(z = l,2,...,c), 



(1) 



(2) 
(3) 



in terms of the energies E and E' , and particle numbers Ni — riiV and = n[V' in both reservoirs. The stationary 
assumption that the pressures, temperatures, and particle densities remain constant in both reservoirs implies that 
the volumes of both reservoirs should satisfy the conditions 



V = - 



and 



V' = + 



(4) 



so that 



T 



V - 



T' 



V' = 0. 



(5) 



Since the system does not exchange entropy with its environment, d^S/dt =^ 0, the time derivative of its total 
entropy S'tot = S + S' gives the entropy production of effusion dStot/dt — dcS/dt + diS/dt — diS/dt. Using Gibbs' 
relation dE — TdS — PdV + t^i^Ni for both gaseous mixtures and Eq. (|5|, the entropy production of effusion 

is given by 



IdiS 
kltt 



= AeJe 



in terms of the affinities or thermodynamic forces defined as 



1 


1 


kT 


~ kT ' 


l-K 




kT 


kT' 



(i = l,2,...,c), 



(6) 

(7) 

(8) 
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where {/Xi}^^]^ and {Mi}i^=i are the chemical potentials of the different species composing the gaseous mixtures in both 
reservoirs. Since the gases are here supposed to be ideal and composed of monoatomic species, the affinities of the 
particle flows are given by 



A, = In 



3/2' 



(i = l,2,...,c). (9) 



If we take the temperature and densities of the reservoir TZ as the reference conditions, the same quantities in the 
other reservoir TZ' can be expressed in terms of the affinities as 

kT' = , (10) 



1 + kTAs ' 
rii exp(-^i) 

(l + /cTA£;)3/2 



(z-l,2,...,c). (11) 



In this regard, the affinities are the control parameters of the nonequilibrium constraints. 

Now, we still need to determine the heat and particle currents ([2])-(|3]). For this purpose, kinetic theory will allow 
us to take into account the mechanistic aspects of effusion. 

B. Kinetic-theory formulation 

Effusion proceeds as a random sequence of particle passages through the pore in both ways between the two 
reservoirs: TZ ^ TZ' . These random events are statistically independent of each other as long as the interaction 
between the particles is supposed to be negligible. In this respect, the size of the pore should be smaller than the 
mean free path of the particles in the gaseous mixture. Therefore, the effects of binary collisions between the particles 
can be ignored and the flow of particles across the pore can be considered to be purely ballistic. We further assume 
that the thickness of the wall separating both reservoirs is smaller than the diameter of the pore so that most particles 
are crossing the pore in straight trajectories with negligible bouncing or interaction inside the pore |31| . 

Under these conditions, the flow of particles can be determined following standard methods from the kinetic theory 
of gases p!8H22] . The mean number of particles passing the cross-section area a of the pore during some time interval 
can be calculated in terms of the Maxwellian velocity distribution and the particle density in the reservoir, from 
which the particles arrive. Since each particle of velocity v carries the kinetic energy e — rriiv'^ /2, the flow of particles 
also determines the heat flow. Kinetic theory describes the transport of energy and particles between the reservoirs 
as a Markovian stochastic process [T8ti22] . The probability Vt{AE, AN) that an energy AE and AN = {ANjYj^^ 
particles are transferred in the direction TZ ^ TZ' during the time interval t is ruled by the master equation: 

-Vt{AE,AN) = V / deWl+\e)Vt{AE~e,A-N~l,) 
at Jo 

+ V / deWt\e)Vt{AE + e,AN + l,) 

POO 

-rt{AE,AN)Y, de Wt\e) + Wt\e) , (12) 
where we use the notation 1^ = {<^jj}j=i with the Kronecker symbol Sij and where the transition rates are given by 

Wt\.) = ^A^li.o-Pl-,^). (14) 



y/2nm,kT' kT' V kT 

for the particle flows TZ ^ TZ' and TZ TZ' respectively [T5H11]. We notice that these transition rates are related to 
each other and to the affinities ([7|-([8]) by 

\ ; ' =exp(e^g + AO (z- l,2,...,c), (15) 
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as the consequence of the nonequihbrium constraints on the flow of particles between both reservoirs. At the equilib- 
rium thermodynamic state where the affinities vanish, the two transition rates are equal and the conditions of detailed 
balancing are recovered. The relations (151 are the analogues for the present system of Schnakenberg's conditions 
defining the macroscopic affinities driving general nonequihbrium stochastic processes [5 ^121 132) . 
We introduce the generating function of the currents and their statistical cumulants as 



Q{Xe,X) = lim --ln(exp(-A£;A£;-A- AN)), , 
where the statistical average is defined by 

(exp(-A£;AS-A- AN))j = ^ Vt{AE, AN) exp {-Xe AE ~ X ■ AN) , 

AB.AN 



(16) 



(17) 



in terms of the probability ruled by the master equation ( [T2| \TT\ [T^ . Taking the time derivative of this average and 
using the master equation, we find in the long-time limit that the generating function is given by 



Q{Xe,X)^Y. W^^+>{e){l-e-^^'-^^)+Wr\e){l-e 



(18) 



Using the transition rates (13) and (14 1, we get the explicit expression of the generating function: 

„-Ai 



QiXE,X) = Y.\'''' 



kT 
2nmi 



1 - 



under the condition that 



(l + kTXE)^ 



1 , 1 



kT' 
27rm, 



1 - 



(1 - kT' Ab)2 



(19) 



(20) 



The expression (19), which was obtained for a gas with a single constituent in Ref. ^22j, is here extended to gaseous 
mixtures. 

The rates of effusion into vacuum are known to be given by 



kT 



I 



= -(Jn,[v, 

ZTrrrii 4 



where 



8kT 



(21) 



(22) 



is the mean value of the speed of particles of mass rrii [TB^ET] . If we use these rates and express the temperature T' 
and the densities {n'^} in terms of the affinities by using Eqs. (lO)-(lll, the generating function can be rewritten as 



g(A£,A) = ^rJl + - 



-Ai 



-(A. -A.) 



(l + kTAE)'^ (l + fcTAs) 



(23) 



[I + kT{AE - Xe)V j 

Consequently, we observe that the generating function satisfies the following symmetry relation: 

QiXE,X) = Q{AE~XE,A~X), (24) 
which is the fluctuation theorem for the currents |S1 • We notice that the generating function vanishes according to 

Q(0,O) = Q(^B,A) =0, (25) 



because of its definition (16), the normalization condition {l)t = 1, and the symmetry relation (24). An equivalent 
expression of the fluctuation theorem for the currents is given by 



Vt{AE,AN) 
Vt{-AE,~AN) 



exp(AE AE + A ■ AN) for t ^ oo . 



(26) 
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In this form, the fluctuation theorem compares the probabiHties of opposite fluctuations in the effusive transport of 
energy and particles across the pore and shows that the ratio of these probabilities is controlled by the affinities. 
At the thermodynamic equilibrium where the affinities are vanishing, the principle of detailed balancing is recovered 
according to which the probabilities of opposite fluctuations are in balance. Out of equilibrium, this balance is broken 
by the noncquilibrium probability distribution and the fluctuating currents tend to follow the directionality fixed by 
the affinities. The fluctuation theorem (26) is the consequence of the relations ( jlsj ) satisfied by the transition rates of 
particles of energy e and species i = l,2,...,c between both reservoirs. The difference between Eqs. (151 and (26 1 is 
that the relations ( 15 1 concern the rates of random transfers of individual particles over an infinitesimal time interval, 
although the fluctuation theorem ( 26 ) holds for the cumulative transport of an energy AE and particle numbers AN 
over the long time interval t. Otherwise, the latter results from the former in much the same way as the fluctuation 
theorem for currents can be deduced from Schnakenberg's definition of macroscopic affinities in general noncquilibrium 
stochastic processes [HHUl [32] ■ 

Now, the average values of the currents can be deduced from the generating function ( 23 ) according to 

dQ 



Je 



dQ 



(0,0) 



lim -{AE)t. 



1 



(0, 0) lim - (AiV, 



(i = l,2,...,c). 



d\i ^ ' ' t-i-oo t 

The energy and particle currents are thus given in terms of the affinities according to 



Je = 2kTj2'' 



1 - 



(1 



-A, 



{1 + kTAEr 



kTAEf 



l,2,...,c). 



(27) 
(28) 

(29) 
(30) 



We recover the well-known result that, if there is the vacuum in the reservoir TZ\ the energy and particle currents 
become 



lim Je 

A— >oo 



2fcr 



4=1 



lim Ji — Ti 

A— s-oo 



(i = l,2,...,c) 



(31) 
(32) 



where is the effusion rate given by Eq. (21 ) 

Using the inequality (e~'^) > e~^'^^ |33j. we have that 



-Ae ab-a-an 



-Ae (AE)t-A-(A-N)t 



(33) 



Therefore, in the macroscopic steady state reached in the long-time limit, we infer from the inequality (33) and 
Eq. (l25l) that 



l<kS_ 
k~dt 



^ AeJe+A-J 
> lim In (e 

t-i-oo t ^ 



lim hAE{AE)t+A-{A-N)t) 

t— yoo t 

-AeAB-A-ANx _ 



)^^Q{Ae,A)^0, 



(34) 



which is the verification that the second law of thermodynamics holds in the present formulation. 



III. NONLINEAR RESPONSE COEFFICIENTS 



The authors of the present paper have previously proved that the symmetry relation ( |24| has for consequence not 
only the Onsager reciprocity relations which only hold in the linear regimes, but also generalizations of these relations 
which extend to the nonlinear response coefficients |11H13[ I36j . These generalizations of Onsager reciprocity relations 
are remarkable, in particular, because they can be used to express the response coefficients in terms of quantities 
characterizing the diffusivities of the currents and higher cumulants in the regime of nonlinear transport by effusion. 
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In this section, we use the following notations: 



A = {Ac,} = {Ae,A] = {Ae, Ai,A2, Ac} : 

A — {Aq} — {A^;, A} = {A^;, Ai, A2, Ac} , 



(35) 
(36) 



with the indices a e {E, 1, 2, c}. 

The response coefficients -/j^./S: M^^p-y^ N^^p-ySj etc... are defined by expanding the average values of the currents in 
powers of the affinities: 



MA) = ^ = La^jjAp + ^Y1 Mc,f3^ApA^ Na^p^sApA^As + ■■■ 



(37) 



Here, the comma separates the index of the current from the indices of the affinities. We notice that the tensors 
Ma,i3'Y, Na^p'yS,... are a priori only symmetric under the exchange of the indices beyond the comma. 

On the other hand, the diffusivities and the higher cumulants of the fiuctuating currents are defined by successive 
derivatives of the generating function with respect to the parameters A: 



Dap{A) E 

Cap-yiA) 



1 d^Q 



2 dXc^dXp 



A=0 



dXadXpdX^ 



A=0 



1 



2 dXadXpdX^dX^ 



(38) 
(39) 
(40) 



so that all these tensors are totally symmetric. 



The fluctuation theorem for the currents ( 24 ) has for consequence the following realtionships starting from the 



Onsager reciprocity relations for the linear response coefficients and extending to the nonlinear response coefficients 

[IIHSIES] 



[ dA, dAp , 



dA^dAs dAadAs dAndA^ 



(41) 
(42) 

(43) 



A=0 



A generalization of Onsager reciprocity relations is given by the total symmetry of the following fourth-order tensor: 



Na,PtS 



\dA^dAs dApdAs dApdA^ 



(44) 



7/ A=0 



which is the consequence of Eq. (43 1 and the total symmetry of the fourth-cumulant tensor (401 [TTI - fT5] . Moreover, 
the fluctuation theorem for the currents ( 24 1 has also for consequence the following relations 
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Bap^l&{^) 



dCg 



OA, 



(45) 



A=0 



Our purpose is now to apply these relationships to the effusion process characterized by the generating function 



(I23f. 

Because the currents are fluctuating, the energy Ai? and the particle numbers AN transferred over a time interval t 
undergo random walks around their mean values, which are characterized by the diffusivities or second-order cumulants 
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(38): 



DEE = i{kTfY.'' 



{l + kTAEY 



Dei = kTri 



(l + kTAE) 



{l + kTAE)\ 



while Dij — for i ^ j. Among the large-deviation properties, we also have the third-order cumulants (39): 



Ceee = 24 (kT)^ J2 ' 



1 - 



[l + kTAEf 



CEEr = 6 {kTf 



1 - 



Ceu = 2 kTri 



l + kTAEY 



{l + kTAEf 



[l + kTAE) 



(46) 

(47) 
(48) 

(49) 

(50) 
(51) 
(52) 



while the other third-order cumulants are vanishing. The fourth-order cumulants can be obtained by using their 
definition ( 40 ) in terms of the generating function or thanks to the relations ( 45 ) : 



dC 



EEE 



dAi 



(0) 



120(A:r)4^r,, 



1=1 



dC 



EEE 



dA, 

QCeEi 

~dA~ 
OCeu 



(0) 



dC 



EEi 



dAf 



(0) =24{kTfn 



(0) = ^(0) ^GikTfn, 
(0) = ^(0) 



dA, ' ' OAe 



^2kTi 



dA, 



-iO) 



(53) 

(54) 
(55) 
(56) 
(57) 



Beeee{0) = 

BEEEi{0) - 
BEEii{0) - 
BEiii{0) - 

Biiii{0) - 

while the other cumulants are vanishing. 

The remarkable relationships (41)-(43) for the response coefficients of the currents (29) and (30) are satisfied as 
well: 



(58) 

(59) 
(60) 
(61) 



Le,e — 


Dee{0) 


^6 {kT f ^r, 

1=1 


LE,i — 


DE^{0) 


^2kTr,, 


Li^E — 


DEriO) 


^2kTr, , 


Li^i — 


A.(0) 


= n, 
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Me, EE 

Ms.Et 
Mi,EE 

Mem 

M,,,E 

M,,u 



^^DE^ 



OAe 
BDe^ 
dA,, 



(0) 
(0) 



^(0) + ^(0) 
OAe^ ' OA, ^ ' 

dD. 



-24(fcT)3 ^r,, 

i=l 

-6 {kTf n , 
-2fcT r,, 
-2 fcT r,; , 



OA, 



(62) 

(63) 
(64) 
(65) 
(66) 
(67) 



E,EEE 



N, 



E,EEi 



N,, 



EEE 



NE,Eii 
Ni,iEE 
NE,m 
Ni,iiE 



^-^{0)-Ibeeee{0) 



^ 2 



EE 



OAeOA, 



(0) + 



Ei 



dAl 



(0) 



1 



B 



EEEi 



= 120 (fcT)-* 



(0) ^24{kTfn 



3Sl^(0) 



dAl 



B 



EEE' 



:(0) 



^(0) 



^(0)-^i?i=;.,.(0) 



= 2 



OAeOA 



(0) 



^(0)-^i?^..(0) 



1 



24 {kTf n , 

&{kTfn, 

6{kTfn, 

2kTn, 

2kTn, 



(0)- -B,,,,(0) 



dAf 2 

which shows that the response coefficients can indeed be obtained from the cumulants. 



(68) 

(69) 
(70) 
(71) 
(72) 
(73) 
(74) 
(75) 



The equahty between the coefficients (591 and (60) is the Onsager reciprocity relation Lea = Li,E, which results 



from the total symmetry of the diffusivity tensor (38 1. However, we point out that, for instance, the equalities between 



the coefficients (63) and (64) [or between the coefficients (73) and ([74])] are specific to the effusion process and they 
are not the consequences of Eqs. (42 )-(45 ) we are here concerned with. A system where such system-specific equalities 



do not hold has been described in Ref!~T4]. Nevertheless, the equalities 



iVB,B£.-2|^^(0)= N,,EEE-2^^^{0) =0, 



E,Ei 



E.iii 



dAEdA; 

cPDee 
dA\ 

, ^''DE^ 

dAj 



(0) 



(0) = N,,uE 



dAl 



dAEdA, 



(0)= 3{kTfn 



(0) -0, 



(76) 
(77) 
(78) 



do result from the total symmetry of the fourth-order tensor ( 44 ) and are the mentioned generalizations of Onsager 



reciprocity relations [TTtiTS] . These equalities, as well as Eqs. (53)- (57), are the consequences of the fluctuation 



theorem (24). We have thus shown that the fluctuation theorem for currents has fundamental implications on the 



nonlinear response coefficients, which indeed satisfy the remarkable relationships (42)- (45) in the effusion process. 
These results hold as well for other effusion processes [23l [24] . 
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IV. MASS SEPARATION 



A. The coupling between the currents 



Here, we consider a binary mixture of particles with different masses toi ^ m2- Mass separation is the process in 
which the molar fraction of one of the species increases as the gaseous mixture flows between the reservoirs. Since the 
mixing entropy should be reduced, mass separation requires energy supply. The effusion process provides this energy 
by the coupling of the particle flows to the heat flow, as shown in previous sections. 

For the effusion of a binary mixture, there exist three coupled flows: the heat flow and the two flows of particles 
of masses toi and m2. In order to achieve mass separation, the difference between the two particle flows must be 
controlled by the other flows. Therefore, we introduce the current of the total number of particles and the difference 
between the particle currents we shall call the differential current |37j 



Jn 
Jd 



■h 

Ji 



J2- 
J2- 



The conjugated affinities are defined as 



An = 



An = 



-(Ai+A2) = -ln 



ni7i2 



1 . , . ^ i , rii /no 



The affinity associated with the differential current is related by Ajj 

^ n[/n'2 
ni/n2 

for the enrichment of the binary mixture by species 1 during the transfer TZ 
The currents are given in terms of the affinities by 

- An -A 

2kTri 



ln(l/v7') to the separation factor: 



Je 



Jn 



Jd 



1 



{i + kTAEr 



2kTr2 



7^'. 



-An+Al 



1 



(l + fcTAg) 



ri 



ri 



(1 



kTAsr 



1 



{1 + kTAE) 



f2 



r2 



^-An+Ai 



1 - 



(1 



-kTAE^ 

-Ajv+Ad 



1 



{i + kTAEr 



If we introduce the molar fractions, also called the molar concentrations. 

Hi 



Til + n2 



(z = l,2) 



(79) 
(80) 



(81) 
(82) 

(83) 

(84) 
(85) 
(86) 

(87) 



the currents can be expressed as follows 

Cl 



Jf 



2a 



C2 



\J2i:m\ \J2'Km2 



PVkT 



\J2'Krai \J2Tim,' 



Jn — o 



Jd = <y 



Cl 






( 


C2 

+ 


\/2'Kmi 


V27rm2y 


Vkf 


\%J2iTmi 


\J2'!im2 


Cl 


C2 ^ 




( 


C'2 


\J2'Km\ 


^/2Trm2j 






\J2Tim2 




(89) 
(90) 



in terms of the pressures and temperatures in both reservoirs. These expressions show that all the currents 
{Je-, Jn-, Jd} are coupled together so that the differential current Jd can be controlled by both the currents Je 
and Jn and, therefore, by the differences of pressures or temperatures. 

The three currents (84)-(86) are depicted in Fig. [2] as a function of the affinity Ae associated with the separation 
factor (83 1. They are compared with their linear approximations, j2™^ = '^B^a,pAp, given in terms of the linear 



response coefficients (58)- (61 1. We observe that the currents manifest important nonlinear behavior and that their 
linear approximation may soon become very crude away from equilibrium. The differential current Jd is observed in 
Fig. [2]: to be positive already in a small interval at negative values of the corresponding affinity Ae- It is in this small 
interval that mass separation occurs because the differential current Jd > is there opposite to the affinity Ad < 0. 
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FIG. 2: The average currents versus the affinity Au in the effusion of particles of masses mi = 1 and m2 — 10 composing 
a binary mixture at the temperature kT = 1 and densities ni — 1 and 712 = 2 through a small pore of area a — 1 from the 
reservoir TZ to the reservoir TZ' at the affinities Ae = 0.1 and An ~ 0.3 with respect to the former: (a) the energy current 
Jb; (b) the total particle current Jm', (c) the differential current Jd- The dashe d lin es depict the linear approximations of the 



currents given by X 



(lin) 



X]fl Lcc^pAp with the linear response coefficients (58l-(61 1. 



B. Fluctuation theorem and entropy production 



Because of microreversibility, the fluctuating transfers of energy AE and particles 



obey the following fluctuation theorem: 
Vt iAE,AN, AD) 



rt{-AE,-AN,-AD) 



AN = Ni+N2., 
AD = Ni~N2, 



' exp{AE AE + AnAN + AdAD) for t ^ 00 , 



(91) 
(92) 



(93) 



and the response coefficients are related to the statistical cumulants by Eqs. (41)-(45). The thermodynamic entropy 



production ( 34 ) is here given by 



1 djS 
k~dt 



^AeJe + AmJn + AdJd>Q, 



(94) 



and is always non negative as the consequence of the fluctuation theorem. We notice that this non-negativity is 
otherwise non trivial to establish, given the analytic expressions of the currents (84)-(86) and the affinities ([t]), (jsTJ, 
and (|82|. 

The positivity of the entropy production under nonequilibrium conditions is illustrated in Fig. [3] where we see that 
it is overestimated by its approximation based on the linearized currents. 

The non-negativity of the entropy production has implications in the organization of the domains where the different 
currents take opposite directions in the three-dimensional space of the affinities. Figure [4] represents two sections in 
this space: the section at A^r = in Fig.J^ and the section at A^; = in Fig. [if^. 

In the plane {A^i^Ae) depicted in Fig. Hli, the directions of the currents Ju and Je switch on the lines Jd = and 
Je = 0. Both lines intersect at the state of thermodynamic equilibrium located at the origin Ad ^ Ae = In the 
vicinity of this point, the linear approximation of the currents holds, but the lines Jd = and Je — become curved 
away from equilibrium because of the nonlinear transport effects. The second law of thermodynamics determines the 
signs of the currents in the different domains limited by the lines Jd = and Je = 0. Indeed, the entropy production 
(94) implies that Ad Jd > on the line = in the plane Ajsi = 0, whereupon we find that Jd > for Ad > 0, as 
it is indeed the case on the right-hand side of the origin in Fig. |4^. As a related consequence, the differential current 
has the opposite sign J^i < on the left-hand side where Ad < 0, again along the line Je = 0- A similar reasoning 
determines the sign of the current Je along the line Jd = 0, as seen in Fig. [4^. Therefore, both currents Jd and Je 
have the same sign in the white domains and they have opposite signs in the grey domains of Fig. [4^, in consistency 
with the second law. 

On the other hand, the entropy production (94) reduces to A^Jn 
Fig. [4jD, where the lines corresponding to Jd 



AdJd > in the plane Ae ^ Q depicted in 
and Jat = are represented. Here also, the different domains are 
organized around the thermodynamic equilibrium state at the origin according to the second law of thermodynamics. 
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FIG. 3: The thermodynamic entropy production ( 94 1 versus the afHinity Ad for effusion in the same conditions as in Fig. [2] The 



dashed hne depicts its approximation given by the hnearized currents ji''"'' — La^pAp with the hnear response coefficients 



(58i-(61i 
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FIG. 4: Effusion of a binary mixture of particles of masses mi = 1 and m2 = 10 at the temperature kT = 1 and densities 
ni = n2 — 1 through a pore of area a — 1: (a) the plane (AdjAe) at Ajv — 0; (b) the plane {Ao, A^) at ^4^ = 0. The domains 
where the currents Jd and Jjv have the same direction are in white and those where the currents have opposite directions are 
in grey 



C. Entropy efficiency 

A concept of entropy efficiency has been introduced by Onsager [251 126) as the efficiency of the thermodiffusive 
separation process to decrease Gibbs' mixing entropy at the expense of increasing the other contributions of entropy. 
In the case of effusion, the entropy efficiency can be defined as minus the ratio between the contribution to entropy 
production caused by the mixing of the two species over the other sources of dissipation: 



AdJd 



AeJe + AnJn 



(95) 



This quantity characterizes the efficiency of mass separation in terms of the power of the differential current divided 
by the power dissipated by the energy and the total particle currents. 

In the regimes where mass separation is effective, the differential current Jd is opposite to the corresponding affinity 
Ao = ln(l/-\/7) so that A^Jd < 0. Therefore, the efficiency (95) is positive because the non-negativity of the entropy 
production (94) imphes that AeJe + A^Jn > —AdJd > 0. Furthermore, the entropy efficiency (95) cannot exceed 



unity, which is its maximum value allowed by the second law of thermodynamics ( 94 1 : 

'7< 1, 



(96) 



in the regimes of effective mass separation where 77 > 0. 
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With regard to the expressions (84)- (86 1 of the currents, we notice that the entropy efficiency (95) depends only on 
four variables: 



T] = ri{kTAE,AN,AD;q), 



(97) 



that are the combination UTAe = (T/T') — 1 of the temperature with the energy affinity, the affinities An and A^, 
and, finally, the ratio of the effusion rates: 



r2 



ni 
n2 



(98) 



which plays a central role. This result shows that the efficiency of mass separation by effusion depends on both the 
mass ratio and the initial concentration ratio. 

In the linear regime close to the thermodynamic equilibrium where the curr ent s can be replaced by their linear 
approximation ji''"^ — La^pAp in terms of the linear response coefficients (58)-(61), the maximum value of the 
entropy efficiency is given by 



/max 



(99) 



which is reached at Ae — and Ad — —Aiq{^—l)/{^Jq + \). Accordingly, the maximum efficiency can approach the 
unit value if either g ^ 1 or g <C 1. However, we should expect deviations from the prediction (99) in the nonlinear 
regimes further away from equilibrium, which is indeed the case. 



0.02 



0.015 



0.01 



0.005 




FIG. 5: The efficiency (95 I versus the affinity An under the same conditions as in Figs. [2]and|3] for which g = 1.58. The dashed 
line depicts the efficiency for the linear approximation of the currents, which would reach the maximum value r/mTx = 0.013 if 
Ae = 0. 



0.02 



0.015 



^ 0.01 



0.005 




FIG. 6: The efficiency (951 versus the affinity Ad for the effusion at the affinities Ae = —0.1 and An = —0.3 of a binary 
mixture of particles of masses mi = 1 and m2 = 10 at the temperature kT = 1 and densities ni = 1 and n2 = 2 in the reservoir 
TZ The pore is of area cr = 1. The dashed line depicts the efficiency for the linear approximation of the currents. Here aiso, 
q = 1.58 and ni^l = 0.013 if Ae = 0. 
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Figure [5] depicts the efficiency ( 95 ) under the same conditions as in Figs. [2] and |3] As seen in Fig. [sj the efficiency 
is positive in a smaU interval of values of the affinity Au extending from the locus where Jd{Au) — to Au = 0. It 
is in this small interval that mass separation is powered by the other currents Je and ■ In Fig. [5] the exact value 
of the efficiency is compared with the prediction of the linear approximation, which is observed to underestimate the 
exact value. However, if negative affinities Ae and Ajv are taken, we observe in Fig.[6]that the linear approximation 
gives an overestimation of the actual efficiency. Since the conditions of Figs. [S] and [6] are close to equilibrium and 
since the ratio ( 98 1 is of the order of unity, the efficiency remains very small with respect to its maximum possible 
value (96). 



0.3 r 



0.2 



0.1 








FIG. 7; The efficiency (951 versus the affinity Ad for mixtures of masses mi = 1 and m2 = 10. The left reservoir TZ is at the 
temperature kT = 1 and densities ni — n2 — 1 while the right reservoir TZ' is at the affinities Ae = 0.5 and A^ = 3. The 
ratio of the effusion rates is here given hy q = 3.16. The dashed line depicts the efficiency for the linear approximation of the 



currents, which would reach the maximum value 77, 



(lin) 



0.078 if Ai 



0. 



Figure [7] shows the efficiency under conditions further away from equilibrium where the effect of nonlinearities is 
stronger. The maximum value of the efficiency is here larger than in the conditions of Fig. [5| Moreover, the actual 
value of the efficiency turns out to be much larger than the value (99) predicted by the linear approximation. We 



remark that the actual value would be smaller than its linear approximation if the affinities Ae and An were negative, 
as already observed in Fig. [6j Accordingly, the largest values of the efficiency are typically obtained for positive rather 
than negative affinities Ae and A]\[. These results clearly demonstrate that the nonlinear dependence of the currents 
on the affinities plays an important role in the process of mass separation by effusion. 




FIG. 8: The entropy efficiency (95 1 for mass separation by effusion versus the affinity An and the mass ratio m2/mi. The left 
reservoir TZ is at the temperature kT — 1 and densities ni = 1 and 712 = 1.5, while the right reservoir TZ' is at the affinities 
Ae = 2 and An ~ 0. The entropy efficiency vanishes at the mass ratio ra^lmx = [n-z/niY ~ 2.25 on the straight line Ad ~ 0. 



The dependence of the efficiency (|95]) on the mass ratio TO2/TO1 and the affinity is depicted in Figs. 8p0 for 
different conditions. In each plot, the efficiency is positive between the straight line Ajj = and the curved line where 
Jd = 0. According to Eq. (86) for the differential current, the intersection between the lines Ao = Q and Jjj = 
happens if the ratio (98) of the effusion rates takes the unit value 9=1 while A^ ^ — 21n(l + kT Ae), in which case 



m2 
mi 



(100) 
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FIG. 9: The entropy efficiency (95l for mass separation by effusion versus the affinity Ad and the logarithm logj^Q(m2/mi) of 
the mass ratio, for the right reservoir TZ' at the affinities Ae ~ 1 and An ~ 10, while the left reservoir TZ is at the temperature 
kT = 1 and densities ni — 1 and 712 — 10. The entropy efficiency vanishes at the mass ratio logj^o('^2/'7ii) = 2 log^Q (722/711) = 2 
on the straight line Ad = 0. Note that the efficiency nearly reaches its maximum possible value (961 allowed by the second law 
under these conditions corresponding to a large difference of pressure between both reservoirs. 




FIG. 10: The entropy efficiency (951 for mass separation by effusion versus the affinity Ad and the logarithm logj^Q(7n2/77ii) of 
the mass ratio, for the right reservoir TZ' at the affinities Ae = 100 and An ~ while the left reservoir TZ is at the temperature 
kT = 1 and densities ni — 1 and 712 — 10. The entropy efficiency vanishes at the mass ratio log (7712/7711) = 2 logj^Q(7i2/7ii) = 2 
on the straight line Ad = 0. Note that the efficiency remains low under these conditions corresponding to a large difference of 
temperature between both reservoirs. 



Figure[8]depicts the efRciency for a moderate temperature difference in a situation where the affinity is vanishing. 
We observe that the efficiency is positive for < if g > 1 and for Ad > ii q < 1. If g = 1 so that the mass 
ratio is related by Eq. ( 100 1 to the initial concentration ratio, mass separation cannot be performed by effusion. Mass 
separation thus becomes effective for g 7^ 1. Nevertheless, the efficiency remains low under the conditions of Fig. [Sj 
which are close to equihbrium. 

Figure [9] shows that the efficiency may approach its maximum possible value ( 96 1 allowed by the second law if the 
affinity An takes large positive values, i.e., for very different pressures in both reservoirs. Indeed, using the expressions 
(84)-(|86| for the currents and the definition (95), the maximum efficiency reached in this limit can be estimated to 



be given by 



q-1 



g + 1 



for A 



N 



(101) 



where q is the ratio defined in Eq. (98). We notice that this maximum efficiency is always larger than the maximum 



value (99 1 of the linear regime. Therefore, the efficiency can be much larger in these nonlinear regimes than predicted 



by the linear approximation. 

In contrast, the efficiency cannot reach high values in the limit where the temperature difference is very large 
because the efficiency decreases as 



Vn 



1 



In(fcTAE) 
kTAR 







for Ae 



-00 . 



(102) 



in this limit, as illustrated in Fig. 
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V. CONCLUSIONS 



In the present paper, we have reported a study of nonhncar transport properties in the process of mass separation by 
effusion in the Ught of the recent advances in nonequihbrium statistical mechanics. Effusion is the baUistic transport of 
particles through a small pore in a wall separating two reservoirs containing gaseous mixtures at different temperatures, 
pressures, and concentrations. This is the classical analogue of ballistic transport in open quantum systems where 
time- reversal symmetry relations known under the name of fluctuation theorem have been recently proved [71 [8] . For 
the effusion of gaseous mixtures, classical kinetic theory based on the Maxwellian velocity distribution can be used to 
obtain the counting statistics and its generating function in terms of the affinities or thermodynamic forces, which are 
the control parameters of the nonequihbrium constraints imposed by the particle reservoirs and driving the transport 
process. 

The fluctuation theorem for the currents established in Refs. [51 [TO] can be considered for the stochastic processes 
of effusion. This fluctuation theorem - which finds its origin in microreversibility - is the consequence of the relations 
(151 between the forward and backward transition rates and the affinities. The main point is that the fluctuation 
theorem for the currents implies that the nonlinear response coefficients obey remarkable relationships that are the 
generalizations of the Onsager reciprocity relations [TTkl3| . as we have here demonstrated for effusion in gaseous 
mixtures. On the one hand, these relationships reveal that the nonlinear response coefficients can be obtained from 
the statistical cumulants characterizing the fluctuations of the currents in nonequihbrium steady states. On the other 
hand, they estabhsh symmetries among the coefficients as given by Eqs. (53)-(57) and Eqs. (76)-(78). 



In section [TV] these considerations are applied to mass separation by effusion in binary mixtures where the flows 
of two particle species are coupled to the energy flow. Here, we have introduced the total particle current and the 
differential current that are coupled not only to the energy current but also together, allowing the differential current 
to be controlled by both the energy and total particle currents. 

The nonlinear dependence of the three currents on the three conjugated affinities is evidenced by comparing their 
actual values to their linear approximations. These results show the importance to go beyond linear response theory 
in considering the nonequihbrium thermodynamics of mass separation by effusion. Because of the nonlinear transport 
effects, the entropy efficiency may be significantly larger than it would be the case in the linear approximation. 

The second law of thermodynamics, which is here proved thanks to the fluctuation theorem, imposes an upper 
bound on the entropy efficiency of mass separation by effusion. This fundamental limit is not reached under moderate 
nonequihbrium conditions, but can nevertheless be approached for large pressure differences between both reservoirs. 

In conclusion, we see that our knowledge of transport phenomena can nowadays be extended from the linear to the 
nonlinear regimes thanks to the recent advances in nonequihbrium statistical mechanics. 



Acknowledgments 

D. Andrieux thanks the F.R.S.-FNRS Belgium for financial support. This research is financially supported by the 
Belgian Federal Government (lAP project "NOSY"). 



[1] L. Onsager, Phys. Rev. 37, 405 (1931). 

[2] D. J. Evans, E. G. D. Cohen, and G. P. Morriss, Phys. Rev. Lett. 71, 2401 (1993). 

[3] G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett. 74, 2694 (1995). 

[4] S. Tasaki and P. Gaspard, J. Stat. Phys. 81, 935 (1995). 

[5] J. Kurchan, J. Phys. A: Math. Gen. 31, 3719 (1998). 

[6] J. L. Lebowitz and H. Spohn, J. Stat. Phys. 95, 333 (1999). 

[7] M. Esposito, U. Harbola, and S. Mukamel, Rev. Mod. Phys. 81, 1665 (2009). 

[8] D. Andrieux, P. Gaspard, T. Monnai, and S. Tasaki, New J. Phys. 11, 043014 (2009). 

[9] D. Andrieux and P. Gaspard, J. Stat. Phys. 127, 107 (2007). 

[10] D. Andrieux, Noneqmlibnum Statistical Thermodynamics at the Nanoscale (VDM Verlag, Saarbriicken, 2009) ISBN 978- 
3639139334. 

[11] D. Andrieux and P. Gaspard, J. Chem. Phys. 121, 6167 (2004). 

[12] D. Andrieux and P. Gaspard, J. Stat. Mech.: Th. Exp. POlOll (2006). 

[13] D. Andrieux and P. Gaspard, J. Stat. Mech.: Th. Exp. P02006 (2007). 

[14] D. Andrieux and P. Gaspard, Phys. Rev. E 77, 031137 (2008). 

[15] R. D. Astumian, Phys. Rev. E 79, 021119 (2009). 

[16] K. Saito and Y. Utsumi, Phys. Rev. B 78, 115429 (2008). 



16 



[17] S. Nakamura, Y. Yamauchi, M. Hashisaka, K. Chida, K. Kobayashi, T. Ono, R. Leturcq, K. Ensslin, K. Saito, Y. Utsumi, 

and A. C. Gossard, Phys. Rev. Lett. 104, 080602 (2010). 
[18] R. Balian, From Microphysics to Macrophysics (Springer, Berlin, 1991). 
[19] M. Knudsen, Ann. d. Physik 28, 999 (1909). 

[20] R. D. Present, Kmettc Theory of Gases (McGraw-Hill, New York, 1958). 

[21] R. K. Pathria, Statistical Mechanics (Pergamon, Oxford, 1972). 

[22] B. Cleuren, C. Van den Broeck, and R. Kawai, Phys. Rev. E 74, 021117 (2006). 

[23] K. Wood, C. Van den Broeck, R. Kawai, and K. Lindenberg, Phys. Rev. E 75, 061116 (2007). 

[24] B. Cleuren, K. Willaert, A. Engel, and C. Van den Broeck, Phys. Rev. E 77, 022103 (2008). 

[25] L. Onsager, Phys. Rev. 55, 1136 (1939). 

[26] R. C. Jones and W. H. Furry, Rev. Mod. Phys. 18, 151 (1946). 

[27] T. De Bonder and P. Van Rysselberghe, Affinity (Stanford University Press, Menlo Park CA, 1936). 
[28] I. Prigogine, Introduction to Thermodynamics of Irreversible Processes (Wiley, New York, 1967). 
[29] G. Nicolis and I. Prigogine, Self- Organization in Nonequihbrium Systems (Wiley, New York, 1977). 
[30] H. B. Callen, Thermodynamics and an Introduction to Thermo statistics (Wiley, New- York, 1985). 
[31] M. Savard, C. Tremblay-Darveau, and G. Gervais, Phys. Rev. Lett. 103, 104502 (2009). 
[32] J. Schnakenberg, Rev. Mod. Phys. 48, 571 (1976). 

[33] This inequality is the consequence of Jensen's inequality {f{X)) > f{{X}) holding for any convex function f{X) and any 
expectation value {•) with respect to some normalized probability distribution such that (1) = 1 [34]. Here, the convex 
function is taken as f{X) — exp(— X). In the context of equilibrium statistical mechanics, this inequality is often called 
the Gibbs-Bogoliubov-Feynman bound [35]. 

[34] L S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic Press, New York, 1980) p. 1101. 

[35] D. Chandler, Introduction to Modern Statistical Mechanics (Oxford University Press, New York, 1987). 

[36] P. Gaspard, in: G. Radons, B. Rumpf, and H. G. Schuster (Editors), Nonlinear Dynamics of Nanosystems (Wiley- VCH, 
Weinheim, 2010) pp. 1-74. 

[37] We notice that we could alternatively introduce the total mass current Jm = mi Ji -|- m2 J2 or the total current Jo = 
+ J2 = P/^/kT — P' /^/kT' , which makes explicit Graham's law. 



